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Abstract. Cosmic rays of kinetic energies below ~1 GeV nucleon"' are thought to play a 
key role in the chemistry and dynamics of the interstellar medium. They are also thought 
to be responsible for nucleosynthesis of the light elements Li, Be, and B. However, very 
little is known about the flux and composition of low-energy cosmic rays since the solar 
modulation effect makes impossible a direct detection of these particles near Earth. We 
first discuss the information that the light elements have brought to cosmic-ray studies. We 
then discuss the prospects for detection of nuclear gamma-ray line emission produced by 
interaction of low-energy cosmic rays with interstellar nuclei. 
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1. Introduction 

Low-energy cosmic rays (LECRs) of kinetic 
energies i 1 GeV nucleon ' are thought to be 
a major player in the process of star forma- 
tion. They are a primary source of ionization 
of heavily shielded, dense molecular clouds 
and the resulting ionization fraction conditions 
both the chemistry in these regions and the 
coupling of the gas with the ambient mag- 
netic field. LECRs also represent an important 
source of heating that contribute to hold molec- 
ular cores in equilibrium against gravitational 
forces. In addition, LECRs are thought to drive 
large-scale magnetohydrodynamic turbulence 
and cause amplification of magnetic field in the 
interstellar medium (ISM; see, e.g., Hanasz et 
al. in these proceedings). 

Despite LECRs being a fundamental com- 
ponent of the Galactic ecosystem, their compo- 
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sition and flux are very uncertain. This is partly 
because these low-energy particles cannot be 
detected near Earth due to the solar modulation 
effect. Although the Voyager 1 and 2 space- 
crafts are approaching the heliopause (i.e. the 
contact discontinuity between the solar wind 
and the ISM) the twin probes will probably not 
be able to observe the local interstellar spec- 
trum of CRs below 1 GeV nucleon because 
of the additional modulation likely occuring in 
the outer heliosheath (Scherer et al. 201 1, and 
Fichtner in these proceedings). 

Thus, for still a long time our knowledge of 
the properties of LECRs will be based mainly 
on analysis of radiation produced by interac- 
tion of these particles with the ISM. In partic- 
ular, a very promising way to study hadronic 
CRs (by far the most abundant in the Galaxy) 
below the kinetic energy threshold for produc- 
tion of neutral pions (280 MeV for p + p colli- 
sions) would be to detect characteristic y-ray 
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Fig. 1. Be abundances vs. [O/H] 
([0/H]=log[(0/H)/(0/H)e], where (0/H)/(0/H)e 
is the O abundance relative to its solar value). Data: 
Rich & Boesgaard (2009); Tan et al. (2009). The 
solid line shows a fit to the data with a secondary Be 
production component plus a primary component 
of constant Be/0 = 1.3 x 10"^ 



lines produced by nuclear collisions of CRs 
with interstellar matter (Sect. 3). However, be- 
fore such a detection will be made possible by 
a significant improvement of the sensitivity of 
the y-ray space instruments, useful information 
on LECRs can already be obtained from mea- 
surements of cosmic abundances of the light 
elements Li, Be, and B (Sect. 2). 



2. Nucleosynthesis of the light 

elements and origin of cosmic rays 

There are two basic reasons to assume that Li, 
Be, and B are significantly produced by nuclear 
interaction of galactic cosmic rays (GCRs) 
with the ISM. First, the abundances of these 
species relative to, e.g.. Si, are much higher 
in the GCRs than in the solar system, by fac- 
tors ranging from 10"* to 10^. Second, as first 
shown by Reeves et al. (1970), the pre-solar 
abundances of these elements as measured in 
meteorites can be approximately reproduced 
by assuming that they are synthesized by frag- 
mentation of interstellar (resp. CR) C, N and 
O by CR (resp. interstellar) protons and a- 
particles, plus a contribution from a -H a re- 
actions for Li production, and further assum- 
ing that the GCR flux and CNO abundances 
have not changed throughout the history of the 



Galaxy'. Quantitatively, Reeves et al. (1970) 
found that a flux of CR protons averaged over 
-10 Gyr Fp{E>3Q MeV)~15 cm'^ s"' would 
be needed to reproduce the pre-solar abun- 
dance of Li (sic). In comparison, the vari- 
ous estimates of the local interstellar spectrum 
(LIS) proposed in the literature over the past 
few decades give Fp(E > 30 MeV) in the range 
~5 - 40 cm"^ s"' (see Indriolo et al. 2009, and 
references therein). 

However, a full calculation of the abun- 
dances of the light elements must take into ac- 
count the CR history and chemical evolution of 
the Galaxy. In such a framework, it was long 
expected that Be and B should be produced 
as secondary species by fragmentation of the 
increasingly abundant CNO nuclei in both the 
GCRs and the ISM, that is that Be and B abun- 
dances measured in stellar atmospheres should 
vary quadratically with the star metallicity. But 
observations made from the 1990's reveal that 
Be and B abundances increase linearly with 
metallicity (or with [O/H], see Fig. 1), at least 
in old halo stars of low metallicity. The ob- 
served B evolution could be accounted for by a 
predominant production of "B by v spallation 
of '^C in core-collapse supernovae (Woosley et 
al. 1990), which is a primary production pro- 
cess. But the observed evolution of Be is still 
not well understood. 



2.1. Nucleosynthesis of primary Be in 
supernova remnants 

Synthesis of primary Be is expected to oc- 
cur in every supernova remnant (SNR) produc- 
ing CRs, because before being released in the 
ISM at the end of the Sedov-Taylor phase, the 
particles accelerated at the blast wave can in- 
teract within the remnant with SN ejecta en- 
riched in freshly synthezised C and O nuclei 
(Parizot & Drury 1999a,b). A simple estimate 
of this Be production can be obtained by as- 
suming that a constant fraction Ocr of the avail- 



' It is now admitted, however, that only *Li, 'Be, 
and '"B are pure products of GCR nucleosynthe- 
sis, 'Li being also significantly produced in the Big 
Bang and in stars and "B by v-induced spallation in 
core-collapse supernovae (see Prantzos 2010). 
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Fig. 2. Evolution as a function of age of the Be pro- 
duction rate per synthesized atom {black curves, 
left axis) and the total luminosity in y-rays > 
100 MeV {red curves, right axis), for a SNR expand- 
ing into a medium of H density ;3ism (lower curves: 
«isM=l cm"^^; upper curves: nisM=100 cm"^). 



able mechanical power processed by the blast 
wave is continuously transformed into kinetic 
energy of CRs, that can produced Be while 
being advected with the plasma downstream 
the shock front. The medium inside the SNR 
is taken to be well-mixed and thus homoge- 
neous. Neglecting in first approximation the 
adiabatic and Coulomb energy losses of the 
fast particles, the Be production rate per O- 
atom in the downstream plasma can then be 
estimated to be {dBe/dt)/0 = ^BeecR, where 
?Be ~ 10"'-' s"' (erg cm""*)"' (O-atom)"' is the 
Be production rate normalized to the CR en- 
ergy density, which post-shock value is ecR - 
kOcRpisuVj. Here, Vs is the forward shock ve- 
locity, pisM the mass density of the ISM sur- 
rounding the SNR and k a factor of order unity 
that depends on the equation of state of the 
shocked gas. We calculated q^e assuming for 
the CR momentum spectrum a phase-space 
distribution resulting from diffusive shock ac- 
celeration, /cr(p) with 4 < s < 4.5, 
and for the relative abundances of the CNO 
elements in the SN ejecta C:N:O=0.2:0.05:l 
(Woosley & Weaver 1995). 

Calculated Be production rates are shown 
in Fig. 2 as a function of time after the SN 
explosion. The evolution of Vs is taken from 



the blast wave solution of Truelove & McKee 
(1999), assuming for the supernova energy, the 
mass of the ejecta, and the power-law index of 
the outer ejecta density profile: /ssn = 1 .5 x 
10^' erg, Mej = 10 Mq and n - 10, respec- 
tively. Also shown in Fig. 2 is the total lu- 
minosity of the SNR in y-rays > 100 MeV, 
which we calculated from the same simple for- 
malism but using the y-ray emissivity qy - 
5 X 10"'"* s"' (erg cm"^)"' (H-atom)"' (Drury 
et al. 1994). We see that both the instantaneous 
Be production rate and the y-ray luminosity in- 
crease with the density of the ambient ISM, 
but also that the blast wave evolution is more 
rapid for larger values of njsM (e g- the end of 
the Sedov phase occurs at frad = 31900 and 
2790 yrs for «ism - 1 and 100 cm"^, respec- 
tively). As a result, the final Be/O abundance 
ratio obtained by integration of the Be produc- 
tion rate over the SNR lifetime obeys the fol- 
lowing scaling relation: 

Thus, a large ambient medium density nisM ~ 
100 cm"^ would be needed to account for the 
observed abundance ratio in low-metallicity 
stars, (Be/0)obs ~ 1.3 x 10"** (see Fig. 1), 
which is not realistic. This robust conclusion 
is in good agreement with the results of the 
more detailed calculations of Parizot & Drury 
(1999a,b). Moreover, as already pointed out 
by Feltzing & Gustafsson (1994), a present- 
day production of Be in SNRs at the level of 
(Be/0)obs would be accompanied by fluxes of 
y-rays > 100 MeV above those observed from 
several SNRs and the Galaxy. 

2.2. On the origin of cosmic rays 

Thus, the observed quasi-linear increase of Be 
abundance with metallicity implies that a sig- 
nificant amount of CNO nuclei in the GCRs 
were - and possibly still are - accelerated out 
of freshly synthesized matter. A first model 
discussed in the literature is that the GCRs 
are accelerated mostly out of SN ejecta rather 
than the average ISM. Thus, Lingenfelter et 
al. (1998) suggested that high-velocity grains 
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formed in SN ejecta could be the source of CR 
metals, which may be accelerated either in the 
forward or reverse shocks of SNe. But the ab- 
sence of radioactive ^'^Ni in the GCRs makes 
this scenario unlikely, because this observation 
suggests a delay i 10^ years between nucle- 
osynthesis of iron-peak nuclei and their ac- 
celeration to CR energies (Wiedenbeck et al. 
1999). Indeed ^^Ni decays by electron capture 
with a half Ufe of 7.6 x 10'' years for the neu- 
tral atoms, but this decay is suppressed once 
the ^^Ni ions are accelerated and fully stripped. 

Another scenario is that the CNO isotopes 
of the GCRs are mainly accelerated from the 
fresh wind material of massive stars, especially 
Wolf-Rayet stars, as it is hit by the blast wave 
of the subsequent SN explosions. Prantzos 
(2010) pointed out that models of fast rotating 
massive stars predict an ejection of winds en- 
riched in CNO elements even at low metalhc- 
ity and suggested that these wind sources could 
produce essentially constant CNO abundances 
in the GCRs at all times. According to Ramaty 
et al. (1997), the measured Be abundances can 
be explained if each Galactic SN supplies on 
average a total energy of ~10^'' erg to GCRs 
of current epoch composition. Assuming a typ- 
ical CR acceleration efficiency of ~10%, it 
means that the SN blast wave has to process 
a large fraction of the total explosion energy 
(~10^' erg) while interacting with a medium 
enriched in freshly synthesized CNO nuclei. In 
other words, the energetics of Be production 
impUes that, on average per Galactic SN, the 
mass of the CNO-rich material ahead of the 
forward shock cannot be much lower than that 
of the SN ejecta. Models for non-rotating mas- 
sive stars do not predict a significant ejection of 
CNO-rich winds at low metallicity. Those for 
fast-rotating, metal-poor, massive stars do dis- 
play a substantial wind enriched in He-burning 
products after the main sequence, but only for 
stars of inital mass > 60 Mq (Decressin et 
al. 2007). For the Salpeter's initial mass func- 
tion, such very massive stars represent no more 
than ~5% of all stars giving rise to SNe. We 
also note that the SN explosion of a fast spin- 
ning star could favor matter ejection through 
jets aligned along the rotation axis (Decressin 
et al. 2007), such that the resulting blast wave 



may not interact with all the wind material. For 
these reasons, it is not clear to us that the sce- 
nario recently put forward by Prantzos (2010) 
can account for enough Be production in the 
early Galaxy. 

A third possibility is that the GCRs are ac- 
celerated primarily out of material inside su- 
perbubbles of hot gas generated by the winds 
and SN explosions of massive stars formed in 
OB associations (Higdon et al. 1998; Parizot 
& Drury 1999c). The majority of core-collapse 
SNe likely occur in these cavities. Shock waves 
propagating in superbubbles can accelerate 
matter enriched in winds and ejecta of SN 
explosions of massive stars born in the same 
OB association, thus producing GCRs of es- 
sentially constant metallicity throughout the 
age of the Galaxy. This model was shown 
to be compatible with the absence of ^^Ni 
in the GCRs (Binns et al. 2008). However, 
the superbubble scenario for the origin of the 
GCRs may suffer from two problems. First, 
contrary to previous results (e.g., Higdon & 
Lingenfelter 2003), Prantzos shows in these 
proceedings that the composition of superbub- 
ble gas can probably not account for the well- 
known anomalous ^^Ne/^^Ne abundance ratio 
in the CR source, which is 5.3 times higher 
than in the solar system. Moreover, specific CR 
acceleration effects are expected in superbub- 
bles due to the presence of multiple shocks and 
plasma turbulence between the shocks (Bykov 
& Fleishman 1992; Parizot et al. 2004; Ferrand 
& Marcowith 2010). As a result, the phase- 
space spectrum of CRs produced in superbub- 
bles is predicted to be harder (resp. softer) 
than /) (the "universal" distribution produced 
by diffuse strong shock acceleration in test- 
particle regime) below (resp. above) a critical 
momentum that depends on the superbubble 
parameters (Ferrand & Marcowith 2010). This 
prediction is not supported by recent Fermi 
LAT y-ray observations (see Mizuno et al. 
2011) that show the spectrum of the y-ray- 
producing CRs to be spatially uniform in the 
outer Galaxy and consistent with the LIS (see 
also Tibaldo et al. in these proceedings). 

It is possible that the bulk of primary Be 
is in fact not synthesized by the GCRs them- 
selves, but by a separate LECR component 
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(Casse et al. 1995) that could be produced 
in superbubbles (Vangioni-Flam et al. 1998). 
These LECRs might be accelerated at early 
times during the superbubble lifetime, when 
the dynamics of the system is powered by 
the strong winds of massive stars. The stan- 
dard GCRs would then be produced later on, 
by regular Fermi (i.e. diffusive shock) accel- 
eration in SN blast waves first propagating in 
the winds of the massive progenitor stars and 
then in the local ISM (in order to explain the 
^^Ne/^"Ne abundance ratio in the GCRs, see 
Prantzos in these proceedings). The postulated 
LECR component could account for the large 
ionization rate deduced from abundances in 
diffuse interstellar clouds (Indriolo et al. 2009), 
although the extra ionization could also be due 
to LECR electrons (Padovani et al. 2009). 
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Fig. 3. Calculated inner Galaxy nuclear y-ray line 
spectrum, with indications of the emitting nuclei for 
some prominent lines. 



3. Nuclear gamma-ray line emission 

The most compelling evidence for a significant 
hadronic LECR component would come from 
the detection of nuclear y-ray fines produced 
by interaction of these fast particles with in- 
terstellar matter. The strongest lines from the 
most abundant nuclei are preferentially pro- 
duced in the CR energy range below a few hun- 
dred MeV per nucleon. We estimate here the 
diffuse Galactic emission in these lines; poten- 
tial y-ray line sources such as SNRs wiU be dis- 
cussed elsewhere. 

Indriolo et al. (2009) proposed several 
LECR spectra which, added to the standard CR 
spectrum, would produce the observed mean 
ionization rate of diffuse molecular clouds 
^CR = 4 X 10"'* s"^ We used a sfightly mod- 
ified version of the broken power-law spec- 
trum of these authors, with the CR composition 
measured at ~1 GeV nucleon"' kept constant 
down to 2 MeV nucleon"', contrary to Indriolo 
et al. (2009) who used solar abundances below 
200 MeV nucleon '. This required a slightly 
changed power-law index of 1 .65 for the LECR 
component to get the same ^cr- It was verified 
that the y-ray emission from decay for such 
an interstellar CR spectrum reproduces the lo- 
cal y-ray emissivity observed by Fermi LAT in 
the range 0.1 - 10 GeV (Abdo et al. 2009). 



For the calculation of the nuclear y-ray 
line emission, cross sections for strong fines 
were taken mainly from Murphy et al. (2009) 
with some modifications taking into account 
new nuclear data of Benhabiles et al. (201 1). 
The production of weaker fines and y-ray emit- 
ting radioactive nuclei was calculated with the 
nuclear reaction code TALYS (Koning et al. 
2008). For the interstellar matter, we used so- 
lar abundances with 2 times solar metaUicity 
to take into account the metal enrichment of 
the inner Galaxy. Parallely to the nuclear lines, 
the y-ray emission following jf' production and 
decay was calculated with the model of Kamae 
et al. (2006). Adjusting the latter to the pion- 
decay spectrum of the inner Galaxy at 0.1 - 
100 GeV as deduced from Fermi LAT observa- 
tions (Strong et al. 20 1 1 ) provided the absolute 
normaUzation. The calculated y-ray line spec- 
trum is shown in Fig. 3. Besides strong narrow 
lines, the nuclear line emission is composed 
of broad lines produced by the CR heavy-ion 
component and of thousands of weaker fines. 
Adding all up results in a complex and ex- 
tremely rich emission spectrum. 

A very promising feature for an observa- 
tion with future space-borne y-ray telescopes 
is the characteristic bump above By ~ J> MeV 
with several strong lines of '^C and '*0. The 
calculated flux in the 3-10 MeV band inte- 
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grated over the inner Galaxy (300° < I < 60°, 
\b\ < 10°) amounts to 7 x lO'^ cm"^ s'K This 
can be compared with the estimated contin- 
uum sensitivities of some recent y-ray tele- 
scope projects, like, e.g., the GRIPS, DUAL 
and CAPSiTT proposals for the latest (2010) 
ESA Cosmic Vision 2015-2025 call for a 
Medium-size space mission. The proposed in- 
struments feature a factor of typically 10 - 
100 sensitivity improvement with respect to 
CG/JO/COMPTEL and INTEGRAL/SPl. Thus, 
the CAPSiTT 3cr continuum sensitivity at 5 
MeV in a 5-year all-sky survey for such a 
spatially extended emission is estimated to be 
~ 10"^ cm"-^ s"\ which would assure a com- 
fortable detection of the predicted nuclear line 
emission. 

4. Conclusions 

The observed primary evolution of Be ver- 
sus [Fe/H] (or [O/H]) shows that some CRs 
are accelerated out of freshly synthesized mat- 
ter. The most promising sites of accelera- 
tion of these CRs appears to be superbubbles, 
which, however, are probably not the source 
of the standard GCRs detected near Earth. 
But superbubbles may produce a distinct com- 
ponent of LECRs during their early, stellar- 
wind-dominated phase. This supplementary 
CR component could account for the large ion- 
ization rate of diffuse interstellar clouds. 

Nuclear interaction y-ray lines could pro- 
vide one of the best ways of studying the origin 
of LECRs and their role in the Galactic ecosys- 
tem. The detection of these lines should be a 
major objective for the next-generation of MeV 
y-ray telescopes. 
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